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FREQUENCY ASYMPTOTES AND ASYMPTOTIC ERROR 
FOR DISTRIBUTED-PARAMETER SYSTEMS 

By W. W. Happ 

NASA, Electronics Research Center 
Cambridge, Mass. 


SUMMARY 

The frequency response techniques for lumped parameter sys- 
tems are reviewed and extended to distributed-parameter circuit 
elements. For example, a three-terminal RC network can be speci- 
fied by a set of asymptotes for amplitude and phase with an 
error usually less than 2 db or 20 degrees respectively. For any 
two— port a preferred set of six asymptotes can be defined from 
which asymptotes for all network function can be derived. This 
set of functions, referred to as deviations, corresponds to net- 
work polynomials in lumped parameter circuits . Properties of 
deviations are siammarized and adapted for computation of frequency 
asymptotes and asymptotic error estimates of distributed-parameter 
systems. Illustrative examples include uniform and tapered thin- 
film networks. 


INTRODUCTION - SCOPE OF INVESTIGATION 
Objectives 

Conventional circuit design concepts and techniques must be 
extended and amended, as distributed parameter networks are re- 
placing lumped parameter networks, in an increasing range of 
circuit applications. Film-type networks not only meet the 
requirements for greater reliability, improved producibility and 
reproducibility, as well as lower cost and often smaller physical 
structures, but in addition have distinct advantages in circuit 
characteristics. For example, distributed R-C structures can 
realize sharper cut-off low-pass filters than lumped R-C circuits. 

As part of design procedures, a circuit designer needs 
analytical expressions representing network-functions for different 
circuit configurations. Distributed parameter networks, such as 
transmission lines, lead to network characteristics described by 
transcendental functions in the complex frequency plane. Appro- 
priate approximations of these transcendental functions provide 
sufficiently accurate design information for thin-film circuits. 
This investigation is directed to investigate the frequency res- 
ponse of distributed-parameter systems by developing: 


(1) An analytical approach to asymptotic approximations 

(2) Design procedures useful for film-type systems 

(3) Universal reference charts or curves for response 
characteristics of thin-film circuit elements. 

Statement of the Problem 

Design techniques for the frequency response in distributed- 
parameter networks are developed by approximating the response of 
the system by a set of asymptotes. This approximation technique 
permits rapid calculation of the response with an error not exceed- 
ing 20 percent. In general, the error introduced by the approxi- 
mation is less than 1 percent. Successive tasks in developing an 
asymptotic description of distributed-parameter systems will be 
as follows: 

(1) Obtain high and low frequency asymptotes for phase and 
amplitude by judicious approximations of transcendental functions 
at the extremes of their range. 

(2) Obtain break-points for asymptotes, i.e., critical 
frequencies separating the regions of high and low frequency 
response. High-frequency asymptotes are valid at a frequency 
exceeding the breakpoint and the low-frequency asymptotes are 
valid at frequencies below the breakpoint. 

(3) Compute order-magnitude estimates of the error intro- 
duced by employing asymptotes for both phase and magnitude. 

(4) Assess the validity of asymptotic approximations by 
a comparison of estimated and computed error. 

(5) Devise a scheme for evaluating the least number of 
functions to describe a three-terminal and eventually four- 
terminal networks . 

(6) Present representative examples illustrating the range 
and effectiveness of method. 

Of the above tasks, (1) and (2) have been successfully completed; 
tasks (4) and (5) have yielded some interesting results; while 
tasks (3) and (6) will be the topics of a future investigation. 

Critical Review of Related Investigation 

Asymptotic techniques are extensively used in lumped para- 
meter networks (refs. 1, 2, 3, 4). To adapt these techniques 
to distributed-parameter networks, it is instructive to review 
the pertinent properties of lumped-parameter asymptotes. The 
network response is described in terms of two singularly inter- 
related components, the amplitude ratio (gain) and the phase- 
difference, with respect to the network response at certain 
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reference frequencies. Approximations of these two functions 
facilitate rapid construction of a graphical presentation suffic- 
iently accurate for engineering calculations. 

Lumped-parameter network functions can generally be expres- 
sed as a ratio of polinomials M(s)/N(s) in n 9 n-negative powers 
of the frequency s = ju, where M(s) = si = Mg tt^ (s-s^) 

= M^; (1 - Tis) and similarly for N(s). Consider now three 

distinct cases: 

(1) If the network has only a finite number R and C (or 

R and L) elements there are, in general, a number of real finite 
roots si = 1/ti = and they yield a set of break frequencies . 
The intersections of frequency asymptotes as obtained on a 
log-log plot are given in Table I. The asymptotes consist of 
straight lines with a discontinuity in the derivative of the 
magnitude and in the phase at the breakpoint. Correction factors 
between actual and asymptotic values for frequencies near the 
break-frequency are given in Table I where the ± depend on oij^. 

On a log-plot both the product of factors (1 - Tis) as well as 
the errors add linearly, permitting rapid compensation of magni- 
tude and phase of M(s) as well as an estimate of the correspond- 
ing error . 

(2) If the network contains R, C, and L elements, complex 
roots Sj^ result in conjugate pairs. Asymptotes and error correc- 
tion curves can be constructed using essentially the same tech- 
nique for expression of the type 1 + as + bs2. 

Table I.- Actual Values and Errors for RC Network 



Range of Validity 

Log Magnitude 

Phase 


0)«^] 


0 

0 

<< 

3 

h— 

< 

=LJ. 

1 


In 

In (w/^i) 

7T/4 

7T/2 

ERROR 

€2 (a>) 


+ 1/2 In (1 + 

1/2 In (1 + -4—,) 

± tan (jj) 

- tan ( _ ) 


3 


(3) If the network contains only energy storing element, 
either L or C or both, roots occur generally at zero of infinity 
and judicious extension of the above techniques are required. 

For film- type circuits only R and C elements are employed, there- 
fore reference is made in Table I to a circuit with an R and a C 
element and unit response at (o = 0. 


REVIEW OF APPLICABLE TWO-PORT THEORY 
Characterization of Frequency Dependent Two-Ports 

Two-ports are characterized by six conventional sets of 
circuit parameters, known as a, b, g, h, y, or z-parameters 
(Table II) . Each of these sets contain four parameters; hence, 
twenty-four parameters need to be considered for a complete in- 
vestigation of a given two-port. Introduction of concept of 
deviations (refs. 5, 6) furnishes a network analyst with the 
following advantages: 

(1) Basic functions characterizing a two-port reduce in 
number from twenty-four to six. For a passive and syinmetrical 
network this number reduces further to four. Hence the concept 
of deviation yields a systematic saving of effort in network 
analysis . 

(2) A circuit designer desiring sufficiently accurate de- 
sign information from asymptotic approximations can easily get 
that by drawing only six asymptotes, which reduce to four asymp- 
totes for passive and symmetrical networks. 

(3) Error analysis in design techniques is similarly reduced 
to a minimum. Because of these advantages associated with devia- 
tions in network analysis, frequency-dependent two-ports will be 
characterized by deviations , Table II shows conventional network 
functions in terms of deviations. To make this presentation self- 
contained and useful to engineers without background in this 
relatively new tool, the rules governing deviations will be brief- 
ly reviewed. Deviations describe multiport networks and are de- 
fined in several alternative ways : 

(a) In terms of cut-sets and tie-sets based on topo- 
logical properties of networks. 

(b) In terms of partial differentials. This approach 
provides a useful definition, but has the advan- 
tage that only ratios of deviations, not deviations 
themselves, are defined. In network applications 
only ratios of deviations are needed, since a par- 
tial derivative can always be expressed as the 
ratio of two deviations. 

(c) In terms of derivatives with respect to a new 
parameter t, which need never be specified. If 
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all the partial derivatives of a system are refer- 
red to some variable t, 


(x,y) 

(u,y) 


If only ratios of deviations are considered, then 
t need not be stated explicitly. This permits 
the definition of the deviation (x,y) as a partial 
derivative with respect to t. 

TABLE II.- NETWORK PARAMETERS IN DEVIATION NOTATION 


Subscri pts 


Deviation 

Symbol 

11 or i 

12 or r 

21 or f 

22 or 0 

Determinant 

(Vq/'q) 

(a) 


_lzL 

( 0 ) 

(°) 

_lbL 

(a) 

_ib} 

(a) 


(b) 

IhL 

(b) 

(b) 

M 

(b) 

(b) 

(a) 

(b) 

(vi,io) 

(g) 

(z) 

(g) 

(b) 

(g) 

(a) 

(g) 

(y) 

(g) 

(h) 

(g) 

('i /^o) 

(h) 

(y) 

(h) 

(b) 

(h) 

(a) 
■ (h) 

(z) 

(h) 

(g) 

(h) 

(vI/Vq) 

(y) 

M 

(y) 

(b) 

(y) 

(a) 

■(y) 

(g) 

(y) 

(z) 

(y) 


(2) 

(g) 

(z) 

(b) 

(z) 

(a) 

(z) 

(h) 

(z) 



9x 

_ 9 X 

/a u 

9u 

y 3 t 

y/* " 


The symbol (x,y) is called the deviation of x with res- 
pect to y. 


To illustrate, conventional h. in Figure 1 



becomes in terms of deviation notation 


h. 


1 




i 


% — ^ 


1 


V 


1 






V 


0 


Figure 1.- Sign conventions on currents and voltages 
in a two-port 

Similarly, other basic rules (ref. 5) have been derived: 

(1) A three-terminal two port is characterized by six devia- 
tions (a), (b) , (g) / (h) , and (z) defined in Table II. 

(2) Table II shows network parameters expressed as ratios 
of deviations. 

(3) Anti-commutation, multiplicative, associative laws 
hold for deviations : 

(x,y) = (-y,x) 

(Kxi,X 2) = K (xi,X2) 

(xi + X2, X3) = (x2,X3) + (X2,X3) 

(4) (xjX2,X3) = XX (X2,X3> + X 2 (xx,X3) 

(5) Any four variables Xx, X 2 f X 3 , X 4 and constrained by 
Xx = f(x 3 , X 4 ) and X 2 = f(x 3 , x^) are related by 

(Xj^,X2) (x^,x^) “H (x2,x^) (Xj|^,x^) t (x2,Xq^) (x2,x^) — 0 

This relationship reduces, for a two-port, to a con- 
straint (a) (b) + (y) - (g) (h) = 0 

(6) An additional deviation (f) is often defined 

(v^ -v^, i^ +i^) = (f) = (g) + (h) - (a) - (b) 

The significance of (f) is seen from the indefinite matrix 
(ref. 7), either [Y] or [Z] , compared in Table III. 
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Evaluation of the elements of the matrices is accomplished 
by grounding a terminal for the [Y] matrix or opening a port for 
the [Z] matrix. For example, Yi2' ^21' ^22 eval- 

uated by grounding the terminal 3. 

Referring to Table II 

yil = h/y y ^2 = 

y2i = -a/y Y22 = g/y 

For convenience of writing deviations , the notation h/y instead 
of (h)/(y) is employed when no ambiguity exists. 

Since the sum of the elements in any row or column of [Z] 
or [Y] is zero, the rest of the elements are found easily. Mat- 
rix relationship can now be written, in terms of deviations as 
given in Table III. 

Distributed Parameter Networks 

The deviations of a simple distributed parameter network 
shown in Table IV are to be derived by proceedings as follows; 

(1) Only deviations (a) , (g) , and (f ) for the common 
capacitor (cc) have to be computed, since for a passive network 
a = b and for a symmetrical network g = h. 

(2) The deviations for the capacitor input (ci) and capaci- 
tor output (co) configurations are obtainable from Table IV using 
the auxiliary variable 

(f) = (g) + (b) - (a) - (b) 

Thus (a) = (b) and (g) = (h) yields 

(g) - (a) = (h)-(b) = (f)/2 

(3) While it is possible to compute deviations directly 

from cut-sets and tie-sets, the deviations (a) , (g) and (f) 

listed in Table III can also be obtained from the matrix develop- 
ment given below. 

(4) For the selection of deviations, an arbitrary constant 
may multiply any deviation. This constant is selected such that 
(a) = (b) = (g) = (h) = 1 at low frequencies for the (cc) network 
in Table III. 
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TABLE 


TWO-PORT 


IV.- DEVIATIONS FOR DISTR 


COMMON CAPACITANCE 
R 

o- wv — 


~r c 


0/SINH 0 
0 COTH 0 


0 COTH 0 
R 

0Vr 


2 0TANH 0/2 


PARAMETER RC NETWORKS 


(CC) I CAPACITANCE OUTPUT (CO)| CAPACITANCE INPUT (Cl) 



0TANH 0/2 

0TANH 0/2 

2 0TANH 0/2 

0 COTH 0 

0 COTH 0 

2 0TANH (0/2) 

R 

R 

2 , 

2 

9/ft 

0 A 

0 COTH 0 

0 COTH 0 






Matrix Development: 




coth0 

-csch0 


(h.) 

-(b) 


e 



1 



Y 

R 



(y) 





-csch0 

coth0 


-(a) 

(g) 


where 0 

R 
C 
A 


/iwCR 

total resistance of the film = r^A 
total capacitance of film = C^A 
length of the film 


This matrix describes the cc configuration and is divided from 
incremental circuit (Figures 2 and 3) such that for 0<x<A 


^ = -ir ^ 

dx o' dx 


Viwc 

o 


^ = Va, i = I where a = 0/A 

dx^ dx^ 


Yielding 


A 

ax 

e 

+ 

B 

e 

ax 

c 

ax 

e 

+ 

D 

e 

ax 


where A, B, C and D are constants to be evaluated from the 
boundary conditions at x = 0 V = I = and at x = A, 

V = V2, I = l2 


Evaluating C and D, 


^1 


0coth 0 

-0/sinh 0 


^1 


1 






R 





^2 


-0/sinh 0 

0coth 0 


^2 


Hence (y) = R; (h) = (9) = 0 coth 0; (la) = (b) = 0/sinh 0 which 
yields (z) = Q^/R. 
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Figure 2.- Distributed Figure 3.- Incremental Model of 

Parameter RC RC Network 

Network 

Asymptotes for Lumped Parameter Network and Deviations 

Frequency domain analysis is a powerful tool in system 
design and is used freguencly in association with the asymptotic 
techniques (refs. 1, 2, 3, 4). 

(1) To determine the stability of a system (Nyquist-Criteria) 

(2) To determine gain - and phase-margin of a given system. 
Alternatively, to design a system with given gain and phase- 
margin. 

(3) To determine the frequency response characteristic of 

an equalizer often used to improve on the response characteristics 
of a given system. 

(4) To determine the close-loop response from a given 
open-loop response with the aid of M and N circles and Nichols 
charts (ref. 2) . 

Frequency response of lumped parameter systems can be des- 
cribed in several alternative ways: 

(1) A polar plot of the frequency response function 

(2) Two complementary curves: one for amplitude of the 

function vs frequency, another for the phase of the function 
vs frequency 

(3) Two complementary curves: one for logarithm of the 

magnitude of the frequency response function vs logarithmic 
frequency scale and another for phase of function vs logarithm 
of the frequency. 

In the discussion to follow, attention is focused on freq- 
uency response. Specifically, properties which can be convenient- 
ly expressed by deviations are: 

(1) The asymptotes may be drawn with a minimum of computa- 
tion . 
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(2) If the true curve is needed, it may be formed by cor- 
recting the asymptotic curve. 

(3) Transfer and driving point functions of lumped para- 
meter networks are rational functions in frequency which can be 
factored. 

Four basic types of frequency functions are used: 


(1) 

K (a 

constant) 

(2) 

(iw) 

± n 

(3) 

(iwT+1) ± n 

(4) 

[ (iw/u^) 2 + 2 


n' 


+ 1 ] 


± n 


Because the logarithm of a function in factored form is the sum 
or logarithms of factors, a composite curve for a frequency res- 
ponse function with many factors can be readily obtained from 
curves corresponding to four basic types mentioned . 


The logarithm of a constant is just another constant and is 
invariant with w. Thus, a plot of 20 log K vs logu is a horizontal 
line. The phase angle is 0 degrees for all frequencies. 

The magnitude curve of the factor (iw)"^ is a straight line 
having a slope of ± 20 n db/decade passing through 0 - db at one 
radian per second. The phase angle is ± nII/2 at all frequencies. 


Magnitude asymptote for the third type of factor consists of 
two parts - a straight line segment coinciding with the 0 db line 
and another- line segment with a slope of ± n db/decade. The 
point of intersection of these Wq asymptotes, known as the break- 
point, is given by coT = 1. Phase asymptotes for this type of 
function are given by two straight lines; one coincides with 
the 0 radian line and another coincides with ± nIl/2 radians line. 
Breakpoints for the phase-asymptotes are made to coincide with 
those of the amplitude asymptotes. 


The fourth type of factor can also be described by asymptotes. 
At low frequencies , the magnitude is unity and the asymptote is 
a horizontal line at 0 db. At high frequencies, the magnitude is 
(w/oon) and the asymptote is a straight line with slope of 
± 20n db/decade. The asymptotes intersect on the 0 db axis at 
w = Wjj, which is the break frequency. Phase asymptotes of this 
type of factor is again represented by two straight lines: one 

coincides with 0-radian axis (for w<Wn) and the other coincides 
with the line passing through ± n (for a)>Wj^) . Magnitude and 
phase asymptotes for different types of factors are given in 
most textbooks (refs. 374). 


12 



It should be noted that factors of the first and second 
types are accurately represented by asymptotes while the third 
and fourth types of factors accrue errors from their asymptotic 
representations . Asymptotic representation of the magnitude of 
third type of function always leads to positive errors (asymptotic 
values greater than the actual values) while errors in the fourth 
type of function depend on the damping ratio § . The true magni- 
tude curve of this type of factor depends on the damping ratio 
which is denoted by §. For §<.707, error is always negative 
while for §>.707, error is positive. Phase asymptotes for the 
third and fourth type of functions always lead to positive errors, 
in the latter case error being a function of § . 

Non-Redundant Network Functions 

The primary motivation for introducing the concept of devia- 
tions is to minimize the number of network functions to completely 
describe a network. For effective presentation of network driving- 
point and transfer functions, the next important step is to 
enumerate non-redundant network functions . In Tables V and VI , 
a three-terminal network is considered and unique driving-point 
and transfer functions, resulting from various configurations 
of the network are listed. 

Driving-point Functions.- To evaluate driving-point functions 
consider the following: 

(1) Exchange of sources (voltage and current) leads to 
reciprocal driving-point functions, impedance and admittance, 
which are not listed separately 

(2) If the network is passive (a) = (b) since the appro- 
priate cut-sets and tie-sets defining deviations are equal for 
passive networks . The reciprocity law in network theory also 
indicates (a) = (b) 

(3) If the network is symmetrical (g) = (h) , since an inter- 
change of terminals, say 1 and 2, transforms (g) into (h) and 
vice versa. This interchange leaves the external characteristics 
invariant and generates no novel network functions. From Table V 
unique driving point functions are: 

(g)/(z) (y)/(h), (f)/(z) and (y)/(f) for the symmetrical 

network , 

and 

(g)/(z), (y)/(h), (h)/(z), (y)/(g), (f)/(z) , (y)/(f), 

for the assymetrical network. 

Transfer Functions.- For the evaluation of non-redundant trans- 
fer functions considerations (2) and (3) in Driving-Point func- 
tions above hold but exchange of sources leads to new transfer 
functions. Additional considerations to be made are: 
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TABLE VII 


DEVIATIONS OF TWO-PORT CONFIGURATIONS RESULTING 
FROM INTERCHANGE OF TERMINALS OR PORTS 
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(1) For a symmetrical network 


'' out in 


= - i . /i . 
out' in 


'out 


out 


(2) For a passive network reciprocity holds. From Table VI 
unique transfer functions are; 

(a)/(g), (a)/(y), (a)/(z), , 


are; 


(h) - (a) 
(h) 


(f) 

for passive, symmetrical networks. 


(3) Transfer functions for passive, unsymmetrical network 


(a)/(g), (a)/(y), (b)/(h), (a)/(z), 

(h) - (a) (h) - (b) (h) - (a) (a) - (g) (a) - (g) 

(y) ' (h) ' ^ ' (g) 

(g) - (b) (a) - (g) 

(f) ' (z) 


Unique driving-point and transfer functions for an active 
network are also obtainable from Tables V and VI assuming 
(a) (b). 

Three-terminal Network.- Transformation of deviations (ref. 8) 
resulting from rotation or reflection of a three-terminal net- 
work are summarized in Table VII using the notation of Table III. 
Note that Table III defines deviations in terms of both, terminal 
and port variables. Terminals and ports transform differently 
under reflection, but the same physical transformation yields 
the same deviations, whether expressed in terminal or port nota- 
tion. Table III in conjunction with Table VII consolidates and 
restates from a somewhat different point-of-view much of the 
information presented in this chapter. 

Representative Distributed Parameter Two-Ports 

Typical values for a, b, g, h, f, y, z for distributed- 
parameter networks are listed in Table VIII for the symmetric 
two-ports and in Table IX for asymmetric two-ports (ref. 9). 

The distributed-parameter RC network is given for comparison with 
Table IV. 


For the transistor (ref. 10) the variables of the two-port 
are taken as currents and minority-carrier concentration densi- 
ties. The resultant parameters are Xp, the minority-carrier 
lifetime, the diff usance 
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H = eADp/W/ and = w^/D^ 

is a time constant defined by w, the width of the base, and D , 
the diffusion constant of minority carriers. The effective ^ 
area is A and the carrier charge is e. 


TABLE VIII.- SYMMETRIC DISTRIBUTED PARAMETER TWO-PORTS 



Thin-f ilm 

Minority-Carrier 

Deviation 

RC Network 

Diffusion Transistor 

a = b 

csch (sRC)^^^ 

csch (ST + T /t ) 
w w p 

11 

coth (sRC)^-^^ 

coth (ST^ + TyXp)^/2 

f 

tanh ^(sRC)^^^ 
1 

tanh |(ST^ + 
1 

y 

(sC/R) 2 

H(1 + sXp)-2 

Z 

(sC/R)^/^ 

(1 + sXp)^/VH 

In Table 

IX the parameter c 

is a measure of the taper, for 


c = 0 the uniform network results. 

TABLE IX.- ASSYMETRIC DISTRIBUTED PARAMETER TWO-PORT 


Deviation 

Tapered RC thin-f ilm network 

a = b 

csch(c^ + sCR)^^^ 

g 

e^coth (c^+sCR) - c/ (c^+sCR) 

h 

e *^coth (c^+sCR) + c/ (c^+sCR)^^^ 

y 

R(c^ + sCR)"^'^^e^ 

z 

r“^(c^ + sCR)^'^^e"‘^ 


FREQUENCY ASYMPTOTES IN DISTRIBUTED-PARAMETER STRUCTURES 

Low and High Frequency Asymptotes for Deviations 

Magnitude and phase asymptotes for the deviations are obtain- 
ed as the limiting values of the function as 0-^0 and 0->«>. As an 
illustration, consider 
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g(0) = 9coth0 


where g(0) is used to denote (g) = (g(0)) consider g(0) as 
0-^0, g(0) = 1 + 0^/3 - 04/45 = ... -»■ 1. Similarly, at 0->-'» 

g (oo) -»-0 

High and low-frequency magnitude asymptotes for deviations 
a and f are derived similarly. Table X gives the asymptotes for 
the three deviations . 

Corresponding phase asymptotes are obtained by evaluating 
phases of the functions at 0-»-O and at 0->-°°. For example, for 
g = 0coth0,,as 0->-O then 

coth0 = 1 
and g = 0 

Similarly for g(0) at 0->-O , the angle Z.g(0) = 0. 

Again, since g(“) = 9 = jojCR, hence 


Z. g = tt/4 


Phase asymptotes for the deviations are shown in Table X 
with (y) = R and (z) = 0/R. 

TABLE X.- ASYMPTOTES FOR DEVIATIONS OF RC NETWORK 


e 

= V^J<JCR = V jx 

Magnitude |f(x)| 

Phase 

^f(x) 

Deviations 

f(x) 

X-3»0 

X— > 00 

x->0 

X— >00 

a = b 

(jx)^/^csch(jx)^/^ 

1 

(4.)'/2„p^(V2) '/2] 

0 

X X 1/2 

V4 - (x/2) 

g = h 

(jx)^/^coth( 

1 

x'/' 

0 

7T/4 

f 

2( jx)^^^tanh(jx/4)^^^ 

X 

! 

2x'/^ 

7T/2 

7T/4 


Note that the deviations listed in Tables IV^ VIII^ and X are only consistent 
within a single table, but otherwise are not unique. 
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Break-Points of Asymptotes 


The intercept of the amplitude asymptotes is termed the 
break-point. The point of intersection is obtained by equating 
the two asymptotes. For example, the intercept x for the asymp- 
totes of a is given by 


/4x 


1 , or X = wRC = 1.05 


Break-points are shown in Table XI. Break-points for phase 
asymptotes are made to coincide with those of magnitude asymptotes 

TABLE XI.- BREAK-POINTS OF DEVIATIONS 


Deviations 

f(0) 


f(0)| 

f(oo) 


Breakpoint 

u 

1 


1 

(2x)l/2(l+i) 

(4x)^/^ 

X = 1 . 05 

a - b 


exp ^x/2{l+i )j 

exp (x/2)V2 

g = h 

1 

1 

{x/2)'/2(l+i) 

(x)’/2 

X = 1 

f 

ix 

X 

(2x)'/^ (Hi) 

2(x)V2 

X = 4 


Magnitude and phase asymptotes tabulated in Table X are plotted 
in Figures 4 and 5. 

Asymptotes of the unique driving-point and transfer functions 
Tablex VIII and IX, can be plotted from Figures 4 and 5 as shown 
in Figures 6 and 7. 

Discussion of Distributed-Parameter Network 
With Tapered Geometry 

R-C distributed parameter networks with rectangular geo- 
metry and uniform resistive and dielectric films have been exten- 
sively investigated and their electrical characteristics are 
generally understood. Attention is now focused on tapered or 
non-uniform d-p networks because of certain desirable electrical 
characteristics not obtainable from rectangular, uniform d-p 
networks. Tapering can be made either geometrically or paramet- 
rically. 
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Heizer (ref. 12) has shown that by proper geometric tapering 
distributed R-C networks with rational transfer functions can be 
realized, thus intensifying the use of these R-C networks (often 
characterized by irrational transfer function) in system design. 
Kaufman (refs. 13, 14) has shown that with exponentially tapered 
RC network null networks can be designed with roll-off charac- 
teristics faster than are obtainable with lumped RC and d-p RC 
circuits. Edson (ref. 15) has proposed tapered RC phase-shift 
oscillators which permit a reduction in the required levels of 
impedance and gain. 

Analyses followed by Kaufman and Edson have been based upon 
parameter tapering only, using a rectangular geometry. 

The case of general tapered RC network, either geometrically 
or parametrically, has not been solved. It has been shown by 
Castro and Happ (ref. 16) that distributed parameter networks 
with circular geometry, which possesses all the desirable charac- 
teristics (Edson, Kaufman) of tapered d-p networks, lend them- 
selves to solution by a one-dimensional formulation. To aid in 
systems design with circular d-p network when applicable, asymp- 
totic plots of short and open-circuit driving point and transfer 
functions and resulting errors are presented. 

Derivation of Deviations For Tapered RC Networks 

An R-C network with circular structure shown in Figure 8 
is equivalent to a tapered network. Although solutions in terms 
of the indefinite admittance and as well as in terms of the in- 
definite impedance matrices, are known, a brief derivation in 
terms of deviations is given. 

Edge effects are not considered; thus, the actual shape of 
the dielectric and conductive films are not considered as long 
as the films extend throughout the resistive film area. The 
outer and inner radii of the resistive film are W and W/n 
(n > 1^ and it is assumed to have a uniform roOhm per square. 

The dielectric film is assumed to have a uniform Cq farads per 
unit area. The matrices are obtained by lumped approximation of 
a distributed circuit as shown in Figure 9. In the Figure 9, dr 
is the elemental resistance of a ring of the film having inner 
and outer radii of 1 and 1 + dl is the elemental capacitance 
associated with it. The voltage and current down the film are 
given by 


^ 

(1) 

d jE 

3^ = iwc^2TTV 

dL o 

(2) 
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noting that 


r 


total 

resistance R 

0 1 

= 27T " 



TT - 

total 

capacitance C 

0 , 2 , , 

2 (n -1) 


and substituting 

R/log n = R', C = 2 n^C/ (n^ -1), 0 = (jwR'C)^'^^ 
Egs . (1) and (2) can be combined to yield 


+ (1/L) (dV/dL) - e^V/W^ = 0 
with the solution 


(3) 


V = AI^ ( 0/WL) + BK^ ( 0/WL) (4) 

where and are modified Bessel functions of first and second 
kind respectively of zero order. Integrating Eg. (2) and apply- 
ing boundary conditions (Figures 2 and 8) 


11 = - |t- ^ [aIj^ ( 0/n) - BK^ ( 0/n)] + D (5) 

1 2 = - { 0/R') [aI^ (0) - BK^(0)] + D (6) 


where and are modified Bessel functions of first order and 
D is an integration constant. Constants A and B are evaluated 
from eguation 


= AI^ ( 0/n) + BK^ ( 0/n) (7) 

V 2 = AI^ (0) + BK^ (0) (8) 
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Define 


A 


1 


I^( e/n) 1^(0) 

K^( 0 /n) K^( 0 ) 


then 

B = [v^I^( 0 ) - V2l^( 0/n)] / Aj_ 


( 9 ) 

( 10 ) 


and constant D equals zero by integrating Eg. ( 2 ) to yield 

V 2 - = f Idr =[a 1^(0) - 0/n)] +[b K^(0) - K^( 0/n)] 


- DR'logn 


( 11 ) 


and substituting Eqs . ( 9 ) and ( 10 ) in Eq. ( 11 ). 

Substituting Eqs. ( 9 ) and ( 10 ) in Eqs. ( 5 ) and (6),. Eqs. ( 5 ) and 
(6) reduce to 


^1 = 

-^1 

A^/A^R' 

- ^2 

Aj/Aj^ 

H 

II 

-^1 

A 4 /Aj^R' 

- ^2 

A5/A1 

where 





^2 = 

0 

Ko(S) 

^0 

( 0/n) 



Io(S) 

^0 

( 0/n) 


I^( 0/n) ( 0/n) 

K^( 0/n) ( 0/n) 


( 12 ) 

( 13 ) 
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^4 = 


^5 = ® 


^0 

(0) 

K'(0) 


(0) 

1^(0) 


( 0/n) 

I' 

0 

(0) 

^0 

( 0/n) 

K' 

0 

(0) 

and (13) 




(y) = R' 
(h) = - • 


(b) 

(a) 

(g) 


A. 


A, 


using the identity 


(x)K^ (x) + (x)I^ 


yields 


"3 = 


and ( a ) = 
hence 


- A, = -1 


(b) = - 


(z) = - 


[It^] 


A^2r. 


(x) = 1 
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Figure 8 . - Incremental model of Circular D-P network 



Figure 9.— A schematic of 
Figure 8 


Asymptotes for Tapered Networks 

The deviations for a circular network, summarized in 
column 1 of Table XII, reveal that a = b and g ^ h, hence, the 
network is asymmetric and passive, as expected. 

Proceedings as in the case of the rectangular network asymp- 
totes at high frequency and low frequency can be found for ampli- 
tude and phase. Proceedings as before, break-points can be calcu- 
lated for specific values of n, or if warranted, universal curves 
for the break-point frequency as a function of n can be plotted. 
Typical values are given in Table XII. The asymptotes are plotted 
in Figures 10 and 11. 


RELATED APPLICATIONS 

To illustrate the usefulness of asymptotic techniques to 
distributed-parameter circuits, the following examples are con- 
sidered: 
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TABLE XII.- DEVIATIONS FOR D-P NETWORK WITH CIRCULAR GEOMETRY 



Magnitude Asymptotes 

Break -Point 

Phase Asymptotes 

Deviations 

f(x) 


)| 

'|x->oo 

arbitrary n 

n=1.5 

^f(x) „ 

X— 


(y) 

-R’ 

-R' 

-R' 



0 

0 

(h) 


l/log n 


nV(log n)^ 

13.7 

0 

+ 74 

(g) 


1/log n 


1 

(log n)^ 

6.1 

0 

+y 4 

(a) 

-^4/A^ 

1/log n 

,„/^,l/2Jx/2)'/^(I-l/n) 


.3 

0 


(b) 

A 3 /A 1 

l/log n 



.3 

0 


( 2 ) 

-(1 +^2 A^)/A^R- 

xK/(logn)^R' 

x/nR' 



V2 

tt /2 

! R‘ = IV^log n 

K = 308/n^ + (1 - 1/n^) .5 log n 






















(1) determination of first null frequencies of four types 
of null filters shown in Table XIII. 

(2) evaluation of transient response of a simple d-p net- 
work (Figure 2 ) . 

Null Frequency.- The first null frequency of each of four filters 
is determined using approximate methods, instead of solving trans- 
cendental Eq. (1). The filter (a) in Table XII is studied in 
detail; results for other networks are tabulated. 

An open circuit voltage-to-voltage transfer function for 
the network (a) is given by 

^21 " '^ 2 /'^! ^1 ^ 

I 2 = 0 

(9/sin h6) + R i (e^R) 

^(0coth0) + ( 0 ^/R)R 2 ^] 

where 

C = Xc^ 

o 

R = Xr^ 

o 

c^ = capacitance per unit length of the d-p network 

= resistance per unit length of the network 
X = length of the film. 

No generality is lost if one assumed that R = C = 1. Null 
frequencies are obtained by setting both the imaginary part and 
the real part of q 2 i simultaneously equal to zero 

Re(g 2 i) and Im(g 2 j^) = 0 

These expressions yield the null frequencies 
tan(w/2)^/^ = -tanh(oj/2)^/^ 

Kaufman (ref. 12) has proposed a graphical solution of this 
equation. 
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TABLE XIII.- NULL FILTERS AND FIRST NULL FREQUENCIES 



921 

Null frequency 
(Radians/sec) 

R] or 

d/sinhe + ( d^/R) 

• 

^cothd +(<?2 /r)r^ 

10.95/(RC) 

R/17.143 

2 

i<jC]d +6 sinh e/R 
j (j C ] 0 cos h 0 + d ^ sinh 0/R 

28.98/(RC) 

6C 

R] 0 + Rsinhd 

0 cos h d + R sinh^ 

28.98/(RC) 

6R 

0/smh0 + i4L> RC] 
B coth 0 + j Ct> RC 1 

10.95/(RC) 

C/17.143 
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Null frequencies, given by points of intersection of 
tan(o3/2)^^^ and -tanh (oj/2) curves, define for which a 

null is obtained from the two equations for 921 ' 

(2)^/^sin(oj/2)^'^^ cosh (to/2) = 1/R^ 

By graphical solution, Kaufman (ref. 14) has found various 
null frequencies and the corresponding values of R^. Our purpose 
is to find the null frequency by approximation technique i.e., 
by approximating the transcendental function to a rational one. 
Only the first null frequency and corresponding value of R^ are 
considered. 

The null is defined by 
9 sinh0 = -R/Rj^ 

By Taylor series expansion 

6(6 + 6^/3! + 0^/5! + 0’^/7l + ) = -R/R^ 


or 

i[co CR-(toCR) Vs'] - I(wCR)V 3I - (toCR)V71] = -R/R^ 
Thus, the first null frequency is given by 
to CR = 10.95 and of course w CR = 0 


and 


R/R^ = (to^CR)V3! - (to^CR)'^/7! = 17.132 

It is interesting to compare these values with those obtain- 
ed by Kaufman (ref. 14) through the graphical approach. His 
results are compared in Table XIV. Thus, this approximation yields 
the desired results with an accuracy within 2.5 percent. Null 
frequencies and the corresponding lumped elements for other fil- 
ters are shown in Table XIII. 
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TABLE XIV.- COMPARISON OF RESULTS 



Graphical 

Asymptotic 

% Discrepancy 

co^CR 

0 

11.18 

10.95 

2.05 

R/Rl 

17.786 

17.143 

3.61 


Transient Response.- Transient response of a simple RC distributed 
network shown in Figure 2 is considered. An open circuit 
voltage-to-voltage transfer function is given 


or 


V2(e)/Vi(0) 



921 


l/cosh0 


V 2 (s)/v^ (s) 


1/cosh (sCR) 


For a step input Ej^(s) = 1/s, so that 


V2(s) 


(1/s) cosh (sCR)^'^^ 


w 1/s [l-sCR/2+5 (sCR) V 24 - 61(sCR)V 720] 
« 1/s - (CR/2 ) [1 - 5sCR/12+61s^C^R^/36o] 

« 1/s -(CR/2) ( ^ - 1+ terms in s^ 


Hence, for t » 1, 

V 2 (t) = l-(6/5,) exp (-12/5CR)t 

At t/CR = 2 

V2(t) « .9897 

This checks satisfactorily with the experimental value .99 given 
by Kaufman and Garrett (ref. 13). A comparison of techniques 
for the transient response is presented by Chen and Bauer (ref. 17) 
and others (refs. 18, 19). 
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CONCLUSION 


The customary presentation of network functions by six con- 
ventional sets of circuit parameters, a, b, g, h, y and 2 has 
been replaced in this investigation by the introduction of the 
concept of deviations. To achieve economy and efficiency in 
evaluating correction curves for a wide variety of distributed- 
parameter network functions, it was necessary to establish a set 
which is sufficient and necessary to derive all possible network 
functions . For a RC two-port it has been necessary to select 
only seven functions referred to as deviations. Evaluation of 
network functions (Tables VII, VIII, and XIII) is accomplished 
by a systematic procedure. Rules for deviations are stated and 
tables presented for reference to accomplish operations of sub- 
stantial complexity. 

Asymptotic techniques for frequency-domain analysis used ex- 
tensively in lumped parameter networks are extended to distributed- 
parameter networks by expansion of deviations at high and low 
frequencies and by establishing "break-points". Delineating the 
range of validity of each asymptotic approximation, frequency 
asymptotes are to be supplemented with error curves to provide 
the engineer with more accurate data for design purposes, in 
case the asymptotic approximations fail to meet his requirements. 

Approximations of transcendental functions by rational func- 
tions give a maximum error of 10 percent as a function of frequen- 
cy. With the use of error-correcting curves an accuracy better 
than 1 percent may be achievable. 

Distributed parameter network functions are transcendental 
functions with an infinite nuinber of ' roots and conversion from 
frequency-domain to time-domain involves evaluation of residues 
at an infinite number of poles. The conventional technique of 
conversion can be replaced by a sufficiently accurate approxima- 
tion obtained by asymptotic expansion of transcendental functions. 
The asymptotic expansion of transcendental network functions 
yield reasonable accuracy in many cases if applied to time-domain 
analysis . 


REFERENCES 

1. Bower, J. L.; and Schultheiss, P. M. : Introduction to the 

Design of Servomechanism. John Wiley and Sons, 1968. 

2. James, H. M.; Nichols, N. B.; and Phillips, R. S.: Theory 

of Servomechanism. McGraw Hill, 1947. 

3. Thaler, G. K. ; and Brown, G. S.; Analysis and Design of 
Feedback Control Systems. McGraw Hill, 1960. 


36 



4. Bode, H. W. ; Network Analysis and Feedback Amplifier De- 
sign. D. Van Nostrand Co., 1945. 

5. Nisbet, T. R.; and Happ, W. W. : Calculus of Deviations 

Applied to Transistor Circuit and Network Analysis. J. Brit. 
Instit. Radio Eng., Vol. 21, 1961, pp. 437-450. 

6. Krajewski, A. S.r and Happ, W. W. ; Deviacje - Nowa Metoda 
Prezeksztatcew Parametrow Elektrycznych, Rosprawy Elektro- 
techneczne, Tom XII, pp. 675-70 6, Warsaw, Sept. 1966. 

7. Shekel, J.: Matrix Representation of Transistor Circuits, 

Proc . Inst. Radio Eng., Vol. 40, Nov. 1952. 

8. Happ, W. W. : Deviations - A Twenty- Year Old Calculus for 

Multitechnical Networks. Proc. Midwest Symposium on Circuit 
Theory, Section IX. 6, Austin, Texas, April 1969. 

9. Happ, W. W. ; and Gupta, S. C^^: Time Domain Analysis and 

Measurement Techniques for RC Structures, J. Appl. Phys . , 

Vol. 40, No. 1, pp. 109-117, Feb. 1969. 

10. Bhattacharyya , B. B.; Stein, R. A.; and Swamy, M. N. S.; 

On Multiwire Non-Uniform Transmission Lines. Proc. Midwest 
Symposium on Circuit Theory, Section II. 1, Austin, Texas, 

April 1969. 

11. Wyndrum, R. W. : Distributed RC Network Synthesis - What 

Important Challenges Remain. Proc. Allerton Conference 
on Circuit and System Theory, p. 286-295, lEE Cat. No. 68 
C44 CT, Oct. 1968. 

12. Heizer, K. W.: Distributed RC Network with Rational Trans- 

fer Function. J.R.E. Transactions on Circuit Theory, 

Dec . 1962. 

13. Kaufman, W. M. ; and Garrett, S. J.: Tapered Distributed 

Filters. I RE Transactions on Circuit Theory , Vol. CT-9 
No. 4, pp. 329-336, Dec. 1962. 

14. Kaufman, W. M. ; Theory of a Monolithic Null Devices and 
Same Novel Circuits. Proc. IRE , Vol. 48, pp. 1540-1545, 

Sept. 1960. 

15. Edson, W. A.: Tapered Distributed RC Lines for Phase-shift 

Oscillators. Proc. IRE , June, 1961. 

16. Castro, P. S.; and Happ, W. W. ; Distributed Parameter Cir- 
cuit Design Techniques. Proc. National Ele ctronics Conference , 
Palo Alto, California, Vol. 17, pp. 47-70, 1961. 

17. Chen, T. H. ; and Bauer, R. F.: On the Transient Response 

of Distributed RC Lines. Proc. Midwest Symposium on Circuit 
Theory, Section XIII, Austin, Texas, April 1969. 


37 


18. Lindholm, F. A.; and Happ, W. W. : Transient Analysis of 

Thin-Film Structures. J. Brit. Instit. Radio Eng., Vol. 26, 
pp. 421-435, 1963. Also Microelectronics and Reliability, 

Vol. 3, pp. 209-226, 1964. 

19. Happ, W. W. ; and Standhammer , J.: Critical Excitation 

Functions of Distributed Parameter Devices. IEEE Trans. 

Vol. ED-13, No. 3, pp. 347-361, March 1968. 

20. Happ, W. W. : Asymptotes for Frequency-Response of Distributed- 

Parameter Components. Proc. Midwest Symposium on Circuit 
Theory, pp. 89-94, Notre Dame, Indiana, May 1968. 


38 


NASA-Langley, 1971 


D-5874 



National Aeronautics and Space Administrai ion 
Washington, D.C 20546 


official business 

PENALTY FOR PRIVATE USE $300 


FIRST CLASS MAIL 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


OIU 001 34 51 3DS 71070 00903 

AIR FORCE WEAPONS LABORATORY /WLOL/ 
KIRTLANO AFBf NEW MEXICO 87117 


ATT E. LOU BOWMANf CHIEFfTECH. LIBRARY 


POSTMASTER; 


If Undeliverable (Section 158 
Postal Manual ) Do Not Return 


^'The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of hutnan knowl- 
edge of phenomena in the atvio sphere and space. The Administration 
shall provide for the xvidest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs^ 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



